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Abstract. The binary parameters of the microquasar GRS 1915+105 have been determined by the detection of Doppler-shifted 
'"CO and '''CO lines in its K-band spectrum (Greiner et al., 2001, Nature, 414, 522). Here, we present further analysis of the 
same K-band VLT spectra and we derive a rotational broadening of the donor star of v sin i = 26 ±3 km s"' from the '"CO/'^^CO 
lines. Assuming that the K-type star is tidally locked to the black hole and is filling its Roche-lobe surface, then the implied 
mass ratio is g = ^ = 0.058 ± 0.033. This result, combined with (P, K, 0=(33.5 d, 140 km s"', 66°) gives a more refined mass 
estimate for the black hole, = 14.0 ± 4.4 Mq, than previously estimated, using an inclination of / = 66° ± 2° (Fender et al. 
1999) as derived from the orientation of the radio jets and a more accurate distance. The mass for the early K-type giant star is 
Md = 0.81 ± 0.53 Mq, consistent with a more evolved stripped-giant donor star in GRS 1915+105 than, for example, the donor 
star of the prototype black-hole X-ray transient, V404 Cyg which has the longest binary period after GRS 1915+105. 
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1. Introduction 

The microquasar GRS 1915+105 remains one of the most ex- 
otic and variable objects in the X-ray, IR and radio spectrum 
bands since its discovery (Castro-Tirado et al. 1994). A major 
accomplishment has been the implied connection between the 
jet and the accretion disk as interpreted through modelling of its 
X-ray variability (Belloni et al. 1997; Klein- Wolt et al. 2002). 
' Thereafter, it is evident that crucial insight into the accretion 
physics of jet formation can be derived through intense theo- 
retical and observational work. Fundamental for the modelling 
but also for the evolutionary history of GRS 1915+105 are the 
binary system parameters. 

The binary model of the microquasar GRS 1915+105 was 
elusive, until recently, despite many observational efforts to 
search for the donor star (Castro-Tirado et al. 1996; Eikenberry 
et al. 1998; Mirabel et al. 1997; Marti et al. 2000; Harlaftis 
et al. 2001; Greiner et al. 2001a). Greiner et al. (2001b) were 
able to obtain phase-resolved IR spectra using the VLT-Antu 
with the ISAAC IR spectrograph and derive the orbital period 
of the binary system from the detected Doppler-shifted '^CO 



and '^CO absorption bands. The deduction of the masses of the 
binary, using the mass function equation. 
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requires the orbital period P, the semi-amplitude of the radial 
velocity curve of the donor star around the black-hole K^, the 
binary inclination / and the mass ratio q - M^/M^, where 
is the mass of the donor star and the mass of the black- 
hole. Indeed, the orbital period and the semi-amplitude of the 
radial velocity curve of the donor star around the black-hole 
have been determined by Greiner et al. (2001b; P - 33.5 + 1.5 
days and Kwa - 140 + 15 km s"'). The remaining unknown 
parameters are the binary inclination and the mass ratio. In 
fact, it is generally the binary inclination which dominates the 
uncertainty in the mass determination of the black hole via 
~ /(Mx) sin ' /. The general method in deriving the inclina- 
tion of X-ray binaries and cataclysmic variables is to model the 
ellipsoidal modulations of the companion star at infrared wave- 
lengths. Luckily, the orbital inclination of GRS 1915 + 105 is 
deduced with unprecedented accuracy from the orientation of 
the jets (/ = 70° + 2°; Mirabel & Rodriguez 1994; / = 66° + 2° 
by Fender et al. 1999). The final parameter that is still undeter- 
mined for GRS 1915+105 is its mass ratio. The mass ratio of 
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the binary system is related to the rotational broadening of the 
donor star using the formula 



y sin; 



= 0.46 [d+qf q] 



1/3 



where v sin / is the rotational broadening in the line-of-sight, 
assuming Roche lobe mass-transfer, spherical geometry for the 
star and tidal locking to the primary star (Gies and Bolton 1986; 
Horne, Wade, Szkody 1986). Here, we further analysed the 
same IR spectra from VLT (Greiner et al. 2001b) in order to de- 
duce the rotational broadening of the donor star from the width 
of the photospheric absorption lines (mainly '^CO, '"^CO; see 
Greiner et al. 2001a). 

2. Observations, numerical technique and 
deduction of mass ratio 

The GRS 1915+105 spectra, together with a Kill spectrum 
(HD 138185), were obtained between April and August 2000 
with the ISAAC spectrograph on the Antu-VLT covering the 
range 2.2896-2.4131 /im, at a pixel resolution of 15.4 km s"' 
(Greiner et al. 2001b). In order to measure the width of the pho- 
tospheric absorption lines, we apply a numerical technique ac- 
cording to which we subtract from the Doppler-corrected spec- 
trum of GRS 1915+105 a Kill template spectrum. In detail 

- we produce the average spectrum after the individual spec- 
tra are Doppler-corrected by applying the radial velocity 
solution (Greiner et al. 2001b) 

- we apply to the Kill template spectrum a rotational profile 
with a limb darkening coefficient of 0.5 (Gray 1992) and 
a veiling factor / which scales the depth of the absorption 
lines in order to simulate the GRS 1915+105 absorption 
lines 

- the difference spectrum represents mainly the accretion 
disk spectrum with all the photospheric absorption bands 
arising from the Kill donor star of GRS 1915+105 having 
been subtracted 

- we perform a x^-test between the difference spectrum and 
its smoothed version (by applying a Gaussian profile of a 
large FWHM=450 km s"' to make it flat) which is a mea- 
sure of how good the absorption lines of the donor star have 
been subtracted for each rotational profile and scaling fac- 
tor </< 1 applied 

- we iterate by changing the rotational profile width and the 
depth of the absorption lines - using a scaling factor / rep- 
resenting the the donor star line flux - until we reach 
values moving along a parabolic surface. 

The minimum on the parabola (v sin ;, x^) gives the 
corresponding width of the rotational profile. We deduce a 
vsini = 26+3 km s"' and / = 14+3 %, at the 90% significance 
level, where / is the light contribution from the donor star as es- 
timated from its '^CO/'-'CO absorption bands. As a validation 
test, we performed the same iterations but now for the individ- 
ual spectra rather than the average Doppler-corrected spectrum 
and took the average v sin / - equal to 26 ± 5 km s ' - from the 
individual estimates. The mass ratio deduced is, only weakly 
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Fig. 1. The spectra show, from bottom to top, the Kill tem- 
plate (with lines multiplied by 0.14), the average Doppler- 
corrected spectrum of GRS I9I5+I05 (16 spectra), the dif- 
ference spectrum (thus, representing the disk spectrum) and 
the difference spectrum between the template and the average 
Doppler-corrected spectrum of GRS 19 15+ 105 at binary phase 
~0.75 (the same average spectrum as in Greiner et al. 2001b). 



sensitive, on the spectral type and veiling factor /. For more de- 
tailed discussion of the technique see Harlaftis et al. (1999 and 
references therein). Fig.[2shows, from bottom to top, the Kill 
template (with lines multiplied by the veiling factor of 0.14), 
the average Doppler-corrected spectrum of GRS 1915+105 (16 
spectra), the difference spectrum (thus, representing the disk 
spectrum) and the difference spectrum between the template 
and the average Doppler-corrected spectrum of GRS 1915+105 
at binary phase ~0.75 (same average as in Greiner et al. 2001b). 
The Kill spectrum has a ripple in the continuum arising from 
variations in the star brightness between integrations (i.e. at dif- 
ferent detector positions). There is a residual blue-shifted CO- 
band or other unidentified emission in the disk spectrum from 
phases other than 0.75 as well as two unidentified absorption 
features around 2.2962 and 2.2988 fim. 

3. Discussion 

We discuss now the various assumptions involved in the deriva- 
tion of the mass ratio from the rotational broadening. The ab- 
sorption vibration-rotation bands '^CO and '^'CO consist of 
~ 4500 lines at 3500 K (Kunde 1977). However the total damp- 
ing width is of the order of 1 x 10"^ A , the thermal width 
does not exceed 0.04 A at 4000 K, and the macroturbulence 
width is ~ 5 km s"' in giant stars (Gray 1992). Other possible 
broadening factors other than rotation is pressure-broadening 
which is not significant in giant stars. Thus, the fact that there 
is a measurable velocity dispersion in the line profiles indi- 
cates that it is dominated by the rotational broadening profile. 
Shahbaz (2003) has developed a new, more accurate method 
involving synthetic spectra from the Roche lobe (thus avoiding 
assumptions about a spherical symmetric star with a limb dark- 
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ening coefficient). However, simulations show that the stan- 
dard method compared to his more accurate treatment is still 
compatible within 3cr. The assumption of a spherical star in 
the (v sin i, q) relation has been shown to introduce an error in 
the evaluation of mass ratio ^ of < 5 % (Marsh et al. 1994) 
which is much less than the uncertainty of our measurement 
(~ 50 %), hence with a neghgible effect. The assumption of 
synchronous rotation (Porbit = ^rot) has been found to be vaUd 
for binaries as wide as 80 days even though giant stars may 
show a weaker synchronism (10% difference from the orbital 
period) than late- type stars (Giuricin et al. 1984). The criterion 
for tidal locking is R^la > 0.125, where is the radius of 
the donor star and a the binary separation, assuming dynam- 
ical tide alone (Zahn 1977; Zahn 1992) which indeed holds 
for GRS 1915-H105 (R^/a = 0.176). Involving other dissipa- 
tive processes such as hydrodynamical meridional flow on the 
star lowers the limit for the fractional star radius, thus yielding 
much longer periods (Tassoul 1987; Tassoul & Tassoul 1992; 
see also for tidal interaction models Witte & Savonije 2002, 
Claret & Cunha 1997, Verbunt & Phinney 1995). The syn- 
chronisation timescale t w (Rd/ay^ q^^ - 4.6 million years 
(Giuricin et al. 1984) which is well within the nuclear evolu- 
tion timescale of the donor star. Synchronisation is more rapid 
than circularisation of the orbit (Hall 1986; see also Lecar et 
al. 1976), and a Xv'test on the radial velocity data (Greiner 
et al. 2001b) does not justify an eccentric orbit (e < 0.006). 
Therefore, we conclude that the orbit is circular. 

The constraint on the eccentricity of the radial velocity fit 
also indicates that there is no Une irradiation affecting the mea- 
surement of the K semi-amplitude of the radial velocity of 
GRS 1915+105. For example, a 10% correction is applied in 
the value of the K semi-amplitude (reduced) from the Nal dou- 
blet absorption lines of the donor star in the eclipsing dwarf 
nova IP Peg from a measured eccentricity of the radial velocity 
fit of (e - 0.089; Martin et al. 1989). There has been no evi- 
dence for irradiation affecting the photospheric absorption lines 
in any quiescent black-hole X-ray transients with extreme mass 
ratios (e.g. Harlaftis et al. 1996). However, outbursts can irra- 
diate the donor star and aff'ect significantly the radial velocity 
amplitude, as for example with the Flll-lV sub-giant donor star 
in GRO J 1655-40 (Shahbaz et al. 2000) which has a binary pe- 
riod of 2.6 days and a large mass ratio of 0.42 (Shahbaz 2003). 
The spectra of GRS 1915+105 were obtained during relatively 
quiet states (no outburst or flare) and in addition its binary orbit 
is long with a small cross section of the donor star for any irra- 
diation to affect it. In fact, using the relation (2) from Phillips 
et al. (1999) we find the effect to be only 2.5% compared to the 
11% uncertainty in Ka for GRS 1915+105. 

The mass ratio thus derived from the rotational broad- 
ening width is q = 0.058 ± 0.033 (see second relation in 
the Introduction). The masses can now be determined from 
iP,K,i,q) =(33.5 ± 1.5 d, 140 ± 15 km s-\ 66° ± 2°, 0.058 ± 
0.033) giving = 14.0 ± 4.4 Mq and Md = 0.81 ± 0.53 Mq 
(see first relation in the Introduction). The large black-hole 
mass uncertainty is dominated by the uncertainty in the semi- 
amplitude radial velocity of 15 km s"^ . The contributions from 
the period uncertainty and mass ratio are about equal but only 
a few % whereas that of the inclination uncertainty is prac- 



tically zero. The large mass ratio uncertainty of 50% is due 
to the large uncertainty in both the rotational broadening and 
the semi-amplitude of the radial velocity (~ 10% for each pa- 
rameter). Having made the above points for the assumptions 
involved we now discuss the ambiguity that may be reside in 
the inclination value. The jet axis may not be perpendicular to 
the binary plane though a constraint is placed to be less than 
25° in GRS 1915+105 (Maccarone 2002). Indeed, in one of 
the microquasars, namely GRO J1655-40, the binary inclina- 
tion is 70° (Greene et al. 2001) whereas the jet inclination is 
found at 85° (Hjellming & Rupen 1995). The proper motions 
of the plasmoids ejected from GRS 1915+105 have been found 
to vary in apparent velocity and this was attributed to a change 
in the jet velocity (Rodriguez & Mirabel 1999; Fender et al. 
1999). However, it must be noted that the alternative explana- 
tion is that the jet axis may be precessing. Indeed, Rodriguez & 
Mirabel (1999) mention a change by ~ 10° of the ejection axis 
within one month as well as a 10-20 % difference in the proper 
motion measured by Fender et al (1999). However, Maccarone 
(2002) concludes, on the basis of close agreement between the 
dynamical estimate and the resonance interpretation of the ob- 
served X-ray QPO frequencies for the black hole mass, that a 
large misalignment between the jet axis and the disk axis is 
ruled out in GRS 1915+105. 

The above discussion must be taken under consideration 
when using the adapted inclination of i = 66° ±2° (Fender et al. 
1999) for a distance of 11 .2 + 0.8 kpc. Assuming a ~ 10° uncer- 
tainty in the inclination due to possible jet axis precession, the 
masses become Mx = 11.6 ± 3.3 Mq and Md = 0.68 ± 0.43 M© 
for an inchnation of / - 76° (M^ = 16.9 + 5.9 Mq and 
Md = 0.98 ± 0.65 Mq for ; = 56°). However, extensive ra- 
dio monitoring of GRS 1915+105 over the last 10 years indi- 
cates no measurable constant precession, thus giving ground to 
the interpretation of varying jet ejection rather than a jet axis 
precession. For direct comparison on the improvement of the 
mass estimate using the mass ratio determined here, we deduce 
a Mx = 12.9 + 4.04 Mq compared to Mx = 14.0 ± 4.0 Mq for 
/ = 70° + 2° as adapted by Greiner et al. (2001a). 

Single Kill stars have negligible usini. Indeed, a KUIl 
and a K5III giant star have a mass of 2.3, 2.2 Mq, a radius 
of 11, 28 ^0 and a < vsini > of 2.5, and less than 1.5 km 
s"', respectively (Gray 1992). Greiner et al. (2001a) identified 
the luminosity class of the donor star in GRS 1915+105 as a 
giant star mainly based on the equivalent width ratio of the 
'^CO/'^CO fines. Using the Eggleton formula (1983) to deter- 
mine the Roche lobe size and hereafter the density of the lobe- 
filUng star, we derive R^ = 19 Rq and p - 2 x 10"^ g cm"^, 
using an orbital period of P = 33.5 days and a mass ratio of 
q = 0.058. This compares well with a single-giant KIIIl-5 clas- 
sification (/? = ll-28/?oandp = 13-4xlO-'*gcm-3).Thus, 
the luminosity class - giant - indicates that the star's envelope 
can be contained within the Roche lobe, and thus evolutionary 
expansion of the donor can sustain Roche lobe mass transfer 
rate. Wind-fed accretion alone cannot power the luminosities 
we observe from this system. 

The long orbital period suggests that the evolution of 
GRS 1915+105 is further advanced than V404 Cyg, a simi- 
lar system with a stripped- subgiant KOIV donor at a period of 
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6.5 days around the black hole (Md = O.Tg jMg; Shahbaz et 
al. 1994). The more evolved giant star in GRS 1915+105 is 
also indicated by the smaller mass of the donor star compared 
to V404 Cyg due to increased mass loss sustained by nuclear 
evolution (King et al. 1996). Indeed, the mass of the donor star, 
a "stripped-giant", is lower by a factor of three compared to a 
single giant star. King (1993) has shown that the mass of the 
donor star can have a maximum mass of I.BMq while its ac- 
tual mass depends on the mass of the helium core. Adapting 
the relation Ra = 12.55/?o(Mc/0.25)5 i from King (1993), we 
derive Mc = 0.27Mo for = 19 Rq. The helium core has a 
mass Mc - 0.21 Mq, one third of the remaining mass of the 
star at - 0.8 IM©. More than 1 Mg mass of the donor 
star has already been transferred to the black hole. The or- 
bital period of GRS 1915+105 is increasing and will reach at 
least 45 days (section 3 in King 1993), a point when the donor 
star will reach its helium-core mass and the mass-transfer rate 
will cease. Actually, the stellar size for the helium-core mass 
of Mc = 0.27Mo is Ri = 18.6 Rq (relation (1) in King 1993), 
suggesting that the steUar size is close to its possible minimum, 
thus not far from the time the mass transfer rate will cease. 

GRS 1915+105 is in perpetual outburst since its discov- 
ery in 1992 whereas V404 Cyg has been in quiescence af- 
ter its X-ray outburst. This is consistent with a higher mass 
transfer rate expected towards the end of the nuclear evolu- 
tion of the donor star. Indeed, a higher mass-transfer rate is ex- 
pected in GRS 1 9 1 5 + 1 05 by 4-8 times compared to V404 Cyg, 
by using the relation, -M2 ^ 5.4 X 10-^(Mc/0.25)'' "Md Mq 
yr"^' (King 1993), compared to an Eddington accretion-rate of 
Mem - 1.5 X 10"^ erg s"^ (for a 14 Mq black-hole mass). 
An increased orbital period, for example from the 6.5 days 
of V404 Cyg to the 33.5 days of GRS 1915+105, due to 
evolution would result in longer outbursts and a decreased 
duty cycle as well as higher super-Eddington accretion rates 
(King et al. 2002). The observed X-ray luminosity during out- 
burst is dependent on inclination and anisotropic radiation pat- 
terns, therefore it is not easy to show that the accretion rate is 
higher in GRS 1915+105. The observed X-ray luminosity is 
super-Eddington at the outburst maximum (Lndd - 1.8 x 10^^ 
erg s"' for a 14 black hole mass) in both V404 Cyg at 
Lx(max) ^ 3.75 x 10^^(0/5 kpc)^ erg s'^ (Tanaka 1989) and 
in GRS 1915+105 at Lx(max) ~ 5 x lO^^iDIM kpc)^ erg s ' 
(Rau et al. 2003). Interestingly, King (2002) recently proposed 
that GRS 1915+105, a low-mass X-ray binary showing super- 
Eddington accretion, may well be the best representative in our 
galaxy of ultraluminous X-ray sources in elUptical galaxies. 

It is clear that, as far as the physical mechanisms involved 
in the radiation emitted, this rich system will be explored more 
efficiently once the binary parameters are known more accu- 
rately. A better sampUng of the orbital period will decrease the 
uncertainty in the stars' masses which mostly arises from the 
radial velocity uncertainty, as well as the independent determi- 
nation of the binary inclination from the IR ellipsoidal modu- 
lations for comparison to the inclination derived from the radio 
jets will constrain further evolutionary considerations of this 
peculiar system. 
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